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Abstract 

Performing lattice simulations of the four dimensional SU{2) gluodynamics we 
find evidence for existence of three-dimensional domains whose total volume scales 
in physical units. Technically, the domains are defined in terms of the minimal 
density of negative links in Z{2) projection of gauge fields. The volume can be 
viewed also as the minimal volume bound by the center vortices. We argue that the 
three-dimensional domains are closely related to confinement. 
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1. It is a general trend in modern theoretical physics to consider extended 
objects, like strings and membranes. Usually, one applies these ideas to hypothet- 
ical, high-dimensional completions of the four-dimensional world. However, lower- 
dimensional structures might also exist in four dimensions. At present time there 
is no well developed theory which would predict such structures. However, there 
is accumulating evidence obtained within the lattice QCD that there are lower di- 
mensions objects percolating through the vacuum of four dimensional Yang-Mills 
theories. We have in mind in particular monopoles and P-vortices, for review see, 
e.g., ^El- The length of the percolating monopole cluster Lp^rc scales in physical 
units PJ: 

Lperc ~ (30.8) (/m)-3 . Vi ; (1) 

where V4 is the volume of the lattice. Similarly, the area of the P-vortices scales in 
the physical units 1^: 

A,ort ~ 24(/m)-2 . F4 . (2) 

Scaling laws , @ have been known since some time but were not orig- 
inally interpreted as an evidence for existence of structures of lower dimensions. 
The reason is that mostly monopoles and vortices were thought of as being 'bulky' 
field configurations with typical sizes of order ^qq£) ■ It is only rather recently that 
it was recognized that at least at presently available lattices they do look actu- 
ally as infinitely thin trajectories and surfaces, i.e. represent physical structures 
of lower dimensions. The basic observation which leads to this conclusion is that 
the monopoles and vortices are distinguished by ultraviolet divergent action, see in 
particular jH EI • In view of this there arises a highly non-trivial question on the 
consistency of the observations with the asymptotic freedom at short distances. The 
data appear to be consistent with the asymptotic freedom 0. 

There is no regular way to search for lower-dimensional defects in the vacuum 
state of a lattice gluodynamics. Historically, the monopoles and P-vortices are 
singled out since they emerged as candidates for confining field configurations. Both 
monopoles and P-vortices are defined in terms of projected fields. In case of the 
monopoles one uses Maximal Abelian projection while in case of P-vortices one 
projects the original SU{2) fields onto the closest Z{2) gauge field configuration. 
The use of a projection makes theoretical analysis on the fundamental level difficult. 
Results ((U are empirical observations. 

In this paper we are looking for possible defects using negative links as the 
scanning means. Namely, we use first a Z{2) projection and then minimize the 
number of the negative links by the residual Z{2) gauge transformations. The 
motivation for such a procedure as well as first encouraging results can be found 
in 0. In brief, negative links correspond to large potentials, A ~ 1/a in the 
continuum limit. Minimizing potentials, on the other hand, might result in gauge 
invariant quantities, for a related discussion see 0. 

We find, indeed, that the negative links, after minimization, occupy a part of 
the lattice which scales as a 3d volume in physical units: 

1/3 ~ 2(/m)-i • F4 . (3) 

To check the projection (in)dependence we perform measurements both in the 
Direct- and Indirect-Maximal Center Projections (DMCP and IMCP). The details 
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of calculations are given in the Appendix. As a result of SU{2) Z{2) projection, 
the original SU{2) field configurations get projected into the closest configuration of 
Z{2) gauge fields. The remaining Z{2) gauge freedom is then fixed by maximizing 
the functional 

F{Z) = Y,Zx,, (4) 



X,fl 



Zx Zx.ii ^: 



±1). In anal- 



with respect to Z(2) gauge transformations {Z. 
ogy with the lattice C/(l) Landau gauge fixing the maximization of the functional 
@ is called ^(2) Landau gauge fixing. 

2. Our results for the minimized density of the negative links are shown in 
Fig. ^ As is seen from the data, both in case of IMCP and in case of DMCP the 
density of negative links is proportional to the lattice spacing. 
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Figure 1: Probability for a link to be negative in Z{2) gauge field configuration after Z{2) 
Landau gauge fixing. 



Let us emphasize that this does imply that there is a strong correlation between 
the negative links and they are accumulating within three-dimensional volumes. 
Indeed, plaquettes dual to P-vortices do have negative links by definition. Thus, if 
we had isolated negative links, then their number would be less than the number of 
negative plaquettes (six negative plaquettes correspond to each isolated link). The 
data indicate, on the other hand, that the number of negative links is larger than 
the number of the negative plaquettes. This could be interpreted only as indication 
that there are many plaquettes with two or four negative links so that they do 
not contribute to the density of the P-vortices but do contribute to density of the 
negative links. The negative link is dual to a 3d elementary cube, the collection 
of these cubes form a 3d volumes which is bound by P-vortices (which also belong 
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to the dual lattice). The Z{2) Landau gauge corresponds to minimization of the 
volumes spanned on Z{2) gauge invariant P-vortices. Thus, it appears that there 
exists a three dimensional volume on the dual lattice whose volume scales in the 
physical units, the fit of the data of Fig. Ogives ©• 

3. Let us mention again that the action associated with the vortices and 
monopoles is ultraviolet divergent jHI^: 

Svor-t ~ 0.53 7, , Smon ~ ln 7 - , (5) 

o 

where S^jort = ~ | < TrC/™''* >) is measured on the plaquettes Up°^* dual to 
P-vortices. The action Smon = ~ ^ < TrC/]?"" >) is measured on the plaquettes 
nearest to monopole currents. From the theoretical point of view it is crucial that 
((Sj) is consistent with the asymptotic freedom [5]. 

However the ultraviolet divergence of the action on the plaquettes belonging to 
the 3d volume would contradict QCD (for a related discussion see 0). We have 
measured this action, S^d, and found that, indeed, 

<Ssd> - <S >raSt"*'^'' ^ , (6) 

in agreement with theoretical expectations. 

4. Now, we will argue that the 3d volumes discussed above might be relevant 
to the mechanism of the confinement. Indeed, let us remind the reader procedure 
introduced in Ref. [H]. One determines Z^^^ variables (see, e.g., the Appendix below) 
and then replaces the original link variables Ux^f^ by new matrices Ux^/j, defined as 

^X,/! = Ux/J, ■ Zx^fJ, , (7) 

where the Z factors are ±1. A remarkable observation of Ref. [Hj is that the string 
tension evaluated by using matrices vanishes and there is no confinement for the 
modified vacuum. 

Usually this procedure is dubbed as 'removal' of the P-vortices. Note, however, 
that the total area of vortices scales in physical units and in this sense they represent 
d=2 defects. Moreover, the Wilson line is obviously a d=l subspace. Thus, in d=4 
the subspaces d=l and d=2 do not intersect at all, generally speaking. For a local 
change (like Q) to affect the Wilson line one should change fields at least on a d=3 
subspace. Thus, no local change of the plaquettes belonging to the vortices can 
eliminate confinement. The resolution of the paradox is that the field modification 
((7)) affects originally links, not plaquettes directly. 

Thus, one can ask what is the minimal number of links which are to be changed 
to eliminate confinement through the procedure ((7|. In this way we come back 
to the 3d volume considered in the bulk of the paper. Moreover, on average the 
number of intersections of a Wilson line with the d=3 volume is given by (see Fig. ^ 
and eq. ©): 

X-'^ inter section) ~ 0.5-t , (8) 

fm 

where Pw is the perimeter of the Wilson line. Note that the number of the inter- 
sections (IHJ does not depend on the lattice spacing a. 
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It is worth emphasizing that if one chooses randomly hnks where the sign of the 
matrix C/^.^^ is changed then the area law cannot be affected and only extra depen- 
dence of the Wilson loop on the perimeter could be generated. Thus, elimination 
of the confinement through multiplication of the Wilson loop by (— l)^™t<!rsection 
possible only in case of a coherent effect. In this context, one can claim these 3d 
defects to be responsible for the confinement. 

Of course the above discussion does not rule interpretation a la Bohm-Aharonov 
effect (which is related to the linking number of the surface with magnetic flux and 
the Wilson line in 4d space ^]). Indeed, it is not possible to eliminate 3d volumes 
without eliminating their boundary, the P-vortices. Moreover the modified Wilson 
loop W = T^Y\i^^{U[} can be expressed through the initial Wilson loop, W, as: 

W' = exp{iTTC{V,C)} W, (9) 

here C is the 4d linking number of the surface of the P-vortex, V, and the Wilson 
contour C. C can be expressed as the number of intersections of C with the 3d 
volumes, or as the integral over V and C. 

Finally let us note that recently 3d manifolds corresponding to the topological 
charge density of the definite sign were found JOl in 4d lattice SU{2) gauge theory. 
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02-04016, INTAS-00-00111, and CRDF award RPI-2364-MO-02. V.I.Z. is partially 
supported by grant INTAS-00-00111 and by DFG grant "From lattice to hadron 
phenomenology" . 



Appendix 

We perform our calculations both in the Direct ^2j - and the Indirect ^Hl ~ Maximal 
Center Projections (DMCP and IMCP). The DMCP in SU(2) lattice gauge theory 
is defined by maximization of the functional 

Fi(C/) = ^(Tr[/,,^)2 , (10) 

X,fl 

with respect to gauge transformations, Ux,f^ is the lattice gauge field. Maximization 
of (fTU]) fixes the gauge up to Z(2) gauge transformations and the corresponding Z(2) 
gauge field is defined as: Z^^^ = signTrC/a,^^. To get IMCP we first fix the maximally 
Abelian gauge maximizing the functional 

F2iU) = ^Tr {U,,,asU+^a,) , (11) 

with respect to gauge transformations. We project gauge degrees of freedom C/(l) — > 
Z(2) by the procedure completely analogous to the DMCP case, that is we maximize 
F\{U) (fTU)) with respect to U{1) gauge transformations. 

To fix the maximally Abelian and direct maximal center gauge we create 20 
randomly gauge transformed copies of the gauge field configuration and apply the 
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Table 1: Parameters of configurations. 





Size 


NiMCP 


Ndmcp 


2.35 


16^ 


20 


20 


2.40 


24^ 


50 


20 


2.45 


24^ 


20 


20 


2.50 


24^ 


50 


20 


2.55 


28^ 


37 


17 


2.60 


28^ 


50 


20 



Simulated Annealing algorithm to fix gauges. We use in calculations that copy which 
correspond to the maximal value of the gauge fixing functional. To fix the indirect 
maximal center gauge from configuration fixed to maximally Abelian gauge and to 
fix the Z(2) degrees of freedom one gauge copy is enough to work with our accuracy. 
We work at various lattice spacings to check the existence of the continuum limit. 
The parameters of our gauge field configurations are listed in Table „ To fix the 
physical scale we use the string tension in lattice units ^Ij, y/a = AAQMeV . 
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